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Abstract 

The e+e" hadron cross section data from PEP, PETRA, TRISTAN, SLC 
and LEP, at centre-of-mass energies between 20 to 209 GeV, are analysed to search 
for the production of a pair of light sbottoms decaying hadronically via R-parity- 
violating couplings. This analysis allows the 95% C.L. exclusion of such a particle 
if its mass is below 7.5 GeV/c^. The light sbottom mass window is closed. 
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1 Introduction 



At the end of the last millenium, the Tevatron Collaborations |2] came out with a 
bottom quark production cross section at ^/s = 1.8 TeV in excess of the theoretical 
prediction by about a factor of two. Refined parton density functions and other theoretical 
improvements, e.g., in the b-quark fragmentation function, have recently been shown to 
account for the difference in the data recorded at = 1.96 TeV 0. 

A more exotic model ||4j, in which a pair of gluinos with mass 12 to 16GeV/c^ is 
produced in pp collisions, with subsequent decays into a bottom quark and a light sbottom, 
with mass below 6GeV/c^, has been shown to also fit the excess well. In this model, the 
sbottom must either be long-lived or decay via R-parity-violating coupling to light quarks, 
e.g., b — > us, to comply with various experimental constraints. Long-lived sbottoms have 
recently been excluded up to masses of 92 GeV/ by ALEPH in direct searches for 
e"'"e~ — > qqqq and e+e~ — >• qq, but R-parity-violating prompt hadronic decays have not 
been addressed by the ALEPH analysis. 

A light, hadronically decaying sbottom would increase the e^e~ — >■ hadron cross section 
above the bb production threshold by up to a quarter of the e"'"e~ bb cross section, i.e., 
about 2% far from the Z peak and 5% at the Z peak. For this reason, the measurements 
of the hadronic cross section at centre-of-mass energies from 20 to 209 GeV {i.e., well 
above the known bb resonances) from PEP PETRA ^ TRISTAN |Ill-|2ni 

and LEP/SLC are re- analysed in this letter to search for a possible consistent excess. 

This letter is organized as follows. A compilation of the data is presented in a synthetic 
manner in Section |21 to allow easy re-interpretation in the future. The global fit of the 
data is described in Section 01 The results of the analysis are given in Section IH and the 
conclusions are listed in Sectional 



2 The hadronic cross section data 

2.1 The data from PEP, PETRA and TRISTAN 

Most of the data from PEP, PETRAN and TRISTAN are published under the form of the 
ratio R of the effective Born hadronic cross section to the point-like e~'~e — >■ n 
cross section 0"° ,, 

_ «Qed(^) ^ 86.85 nb 

where s is the e"'"e~ centre-of-mass energy squared and Oqed is the fine structure constant. 
The latest TOPAZ jTHl Ej and VENUS 20.^ publications report directly the value of a^^^^ 
instead. In both cases, the latter includes a correction that unfolds the effects of initial 
state radiation (ISR), while still reflecting the running of the fine structure constant with 
the centre-of-mass energy [22] • 

The R and a^^^ data are hsted in Tabled (PEP, PETRA) and in Table H (TRISTAN), 
as obtained from a comparison of two recent compilations [23 123 and the original pub- 
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Table 1: The ratio R and the effective Born hadronic cross section, crfJaj, from the PEP and PETRA 
experiments, with increasing centre-of-mass energy (y/s). The expected statistical (astat), point-to-point 
systematic (cptp) and normaUzation systematic (Anorm) uncertainties are also given (in %). The latter is 
correlated between all energy points in a given publication. An additional normalization error Aqed = 
0.1%, fully correlated between all measurements, is to be added to account for missing QED higher orders. 
The last column points to the original publication. 
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Table 2: The ratio R and the effective Born hadronic cross section, <J^g_^, from the TRISTAN ex- 
periments, with increasing centre-of-mass energy (\/s)- The expected statistical ((istat), point-to-point 
systematic (cptp) and normalization systematic (Anorm) uncertainties are also given (in %). The lat- 
ter is correlated between all energy points in a given publication. An additional normalization error 
Aqed = 0.1%, fully correlated between all measurements, is to be added to account for missing QED 
higher orders. The last column points to the original publication. 
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lications 0^1201 • In these tables, only the final - and most accurate - result for each 
experiment and each centre-of-mass energy is reported. (Superseded data are reported in 
both Refs. and |21j, but are not always clearly fiagged as such therein.) 

Other refinements were considered in this letter for a rigorous statistical treatment of 
the data, and are described in the following. First, in each experiment, the systematic 
uncertainty was divided into a point-to-point contribution, (jptp, and an overall normal- 
ization error, Anorm, as is done in most of the orginal publications. The point-to-point 
systematic uncertainties are uncorrelated (related to, e.g., the limited simulated statistics, 
or the statistical uncertainty on the measured luminosity), are assumed to have a Gaussian 
probability density function and are taken directly from the original publications. 

In contrast, the overall normalization error definition varies among the publications, 
being either the largest possible variation interval {e.g., between several sets of selection 
criteria, different ways of determining the luminosity, or various quark fragmentation 
models) or half this interval. Here, the definition was unified in such a way that the 
overall normalization can vary by —Anorm and +Anorm, with a uniform probability over the 
whole interval. This overall normalization error is 100% correlated between the different 
centre-of-mass energy points reported in each given publication. 

Third, the published values of Anorm often contain an estimate of the effect of missing 
higher-order QED corrections in the ISR unfolding procedure, at the level of a couple of 
percent. Indeed, at the time of PEP, PETRA and TRISTAN, the Monte Carlo programs 
used to simulate the e"'"e~ — *• qq and e+e" —>■ e~^e~ processes were limited to 0(q;qed)- 
The missing orders have a potential effect on the measured value of a^^^ via the prediction 
of both the hadronic cross section and the Bhabha scattering cross section: the former 
is used to correct the measured o"had for QED effects, and the latter to determine the 
integrated luminosity. Altogether, the published cross section values would have to be 
corrected as follows, 

„Q ^ ^0 V "cc "had /r,\ 

where the indices (1) and (all) refer to the cross section prediction up to the QED first 
order (used in the original publications) and with all orders, respectively. 

With the programs that have been developed for LEP, it is now possible to evalu- 
ate this correction with a better accuracy than that assumed twenty years ago. The 
e"'"e~ — > qq and e^e~ — >■ e"'"e~ cross sections were determined here with and without QED 
higher orders by ZFITTER [23 with an emulation of the kinematical cuts described in the 
original publications. It was found that the corrections to Bhabha scattering and hadron 
production essentially cancel in the ratio of Eq. 2. The remaining contribution of QED 
higher orders is at the 0.1% level, almost independently of the event selection and the 
centre-of-mass energy. 

The large uncertainties related to the missing QED higher orders were therefore taken 
out from the original values of Anorm- While the aforementioned 0.1% contribution could 
be simply corrected for in cr^^d' ^ normalization error Aqed = 0.1% was added instead 
(assumed to be 100% correlated between all PEP, PETRA and TRISTAN measurements) 
to conservatively account for the yet missing orders in ZFITTER. 
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Finally, early TRISTAN data jTH [T3 dl Ull are also corrected in the original pub- 
lications for other electroweak effects, dominated by the top quark contribution (with 
a {jntop/mzY dependence at first order). These small corrections (between +0.1% and 
+0.7% at ^/s = 60GeV, depending on the top quark mass chosen to determine the cor- 
rection) were unfolded here (i) to have a consistent data set to work with; and (ii) for 
a sound comparison with the ZFITTER prediction, which includes first- and higher-order 
electroweak contributions as well. The latest TRISTAN data [1111131201 were, more ade- 
quately, corrected for QED effects only. The electroweak effect correction needs therefore 
not be unfolded in that case. 

For practical reasons, the measurements of Tables ^ and |21 were clustered in few centre- 
of-mass bins as indicated by the horizontal separation lines in these two tables. The 
ratio R values were averaged in each bin according to the total uncertainties, i.e., with 
a weight proportional to the inverse of a^^^ = x (a^^^ + a^^p + Al^^jS + A^^d/s)- 
The corresponding averaged Born effective cross sections ((7°^^^) and centre-of-mass energy 
values are displayed in Table El The R values found for PEP/PETRA were found to agree 
with those of an earlier combination |H1. The effective Born hadronic cross section (cr^h) 
predicted by ZFITTER [221 is also shown in Table El 

Table 3: The ratio R and the effective Born hadronic cross section, cr{J^j, from PEP, PETRA, TRISTAN, 
as a function of the centre-of-mass energy {^/s), averaged in ~ 2-GeV-wide centre-of-mass-energy bins. 
The hadronic cross-section prediction, cr^jj, is also shown. The last column displays the number of 
measurements used for each entry. 



(GeV) 


Ratio R 


f^had (pb) 


< (pb) 


iVpts 


21.995 


3.843 + 0.067 


754.1 + 13.1 


763.1 
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616.4 + 25.4 
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28.932 


3.945 + 0.045 


450.3 + 5.2 
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30.570 


3.929 + 0.086 


402.2 + 8.8 
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34.408 


3.996 + 0.038 


323.9 + 3.1 
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3.871 + 0.102 


286.6 + 7.5 


291.0 
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38.237 


3.894 + 0.105 


255.9 + 6.9 


260.2 
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41.329 


4.083 + 0.081 


230.3 + 4.6 


225.7 
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43.825 


4.027 + 0.051 


202.2 + 2.6 
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46.038 


4.234 + 0.098 


192.9 + 4.5 
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53.097 


4.527 + 0.097 


155.8 + 3.3 


153.6 


12 


56.432 


4.964 + 0.087 


151.1 + 2.6 


145.4 


9 


57.867 


4.926 + 0.046 


142.7+1.3 


143.4 


16 


60.264 


5.456 + 0.107 


145.8 + 2.9 


142.2 


9 


61.521 


5.378 + 0.156 


138.0 + 4.0 


142.8 
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Computer-readable files for these data will be transmitted to the Review of Particle 
Physics and are available at jhttp : // j anot . web . cern . ch/ j anot/HadronicData/ j 
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The ratio and the difference of these measured cross sections and those predicted 
by ZFITTER are displayed in Fig. ^ as a function of the centre-of-mass energy. When 
no systematic uncertainties are assigned to the theoretical prediction, the average ratio 
appears to exceed the prediction by (0.79 ± 0.52)%, i.e., by 1.5 standard deviations. This 
excess is, however, about 2.4 standard deviations below the prediction of an additional 
light sbottom pair production (here with a mass of 6GeV/c^), which would amount to 
about 2% of the total cross section. 
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Figure 1: Ratio (a) and difference (b) of the effective Born hadronic cross section measurements and 
the ZFITTER prediction as a function of the centre-of-mass energy, for PEP, PETRA and TRISTAN data, 
rebinned as explained in the text. The dash-dotted line indicates the standard model prediction, and the 
dashed curve the additional contribution of sbottom pair production with = 6 GeV/c^ and with a 
vanishing coupling to the Z. 

The experimental correlations between the different bins, essential for a rigourous sta- 
tistical treatment of the data, were determined following the lines of Ref. jHl- In practice, 
a Monte Carlo technique relying on the generation of many gedanken experiments was 
used to determine the probabihty density functions of the measured R ratio values listed 
in Tables ^ and 121 In each gedanken experiment, 108 R values were generated around the 
measured central value, smeared by (i) a, Gaussian distribution with a width equal to the 
quadratic sum of Ustat and cxptp; (ii) a uniform distribution in the [— Anorm, +Anorm] inter- 
val, identical for all energy points of a given publication; and (Hi) a uniform distribution 
in the [— Aqed, +Aqed] interval, identical for all 108 measurements. 
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As above, an average value Ri was determined in each centre-of-mass-energy bin i for 
each gedanken experiment. This allowed the Ri values of Table El and their uncertainties 
to be confirmed, when averaging over a large number of gedanken experiments. Similarly, 
the uncertainties of the cross-products Ri x Rj led to the correlation matrices shown in 
Tables |3] and El for PEP and PETRA on the one hand, and for TRISTAN on the other. 
The cross-correlations between PEP, PETRA and TRISTAN (induced solely by Aqed) 
were found to be smaller than 5 lO""^ and were therefore neglected in the following. 



Table 4: The correlations between the ten PEP and PETRA centre-of-mass energy bins (-^s in GeV). 





21.994 


25.003 


28.932 


30.572 


34.409 


36.027 


38.231 


41.325 


43.824 


46.042 


21.994 


1.000 


0.034 


0.003 


0.043 


0.237 


0.060 


0.028 


0.084 


0.227 


0.066 


25.003 


0.034 


1.000 


0.004 


0.053 


0.096 


0.055 


0.012 


0.053 


0.029 


0.041 


28.932 


0.003 


0.004 


1.000 


0.032 


0.017 


0.009 


0.000 


0.003 


0.001 


0.002 


30.572 


0.043 


0.053 


0.032 


1.000 


0.198 


0.105 


0.009 


0.057 


0.030 


0.043 


34.409 


0.237 


0.096 


0.017 


0.198 


1.000 


0.185 


0.032 


0.148 


0.258 


0.113 


36.027 


0.060 


0.055 


0.009 


0.105 


0.185 


1.000 


0.010 


0.096 


0.048 


0.071 


38.231 


0.028 


0.012 


0.000 


0.009 


0.032 


0.010 


1.000 


0.030 


0.054 


0.030 


41.325 


0.084 


0.053 


0.003 


0.057 


0.148 


0.096 


0.030 


1.000 


0.078 


0.098 


43.824 


0.227 


0.029 


0.001 


0.030 


0.258 


0.048 


0.054 


0.078 


1.000 


0.066 


46.042 


0.066 


0.041 


0.002 


0.043 


0.113 


0.071 


0.030 


0.098 


0.066 


1.000 



Table 5: The correlations between the five TRISTAN centre-of-mass energy bins (-^/s in GeV). 





53.141 


56.436 


57.863 


60.253 


61.519 


53.141 


1.000 


0.101 


0.014 


0.077 


0.057 


56.436 


0.101 


1.000 


0.017 


0.093 


0.073 


57.863 


0.014 


0.017 


1.000 


0.022 


0.010 


60.253 


0.077 


0.093 


0.022 


1.000 


0.058 


61.519 


0.057 


0.073 


0.010 


0.058 


1.000 



2.2 The LEP 1 and SLC data 



Table 6: Precise LEP and SLC measurements of the Z lineshape parameters (Fz, Ri, (Jhad), of gv/dA 
and of i?b, together with their correlation matrix. The last two measurements have been taken here as 
uncorrelated with the first three ^ . The standard model prediction formula are given in Ref. |27| . 



Observable 


Measurement 


Correlation matrix 


Tz 


2495.2 ± 2.3 MeV 


1.000 


Ri 


20.767 ± 0.025 


+0.004 1.000 


Chad 


41.541 ± 0.037nb 


-0.297 -f0.183 1.000 


9v/gA 


0.07408 ± 0.00068 


0.000 0.000 0.000 1.000 


Rh 


0.21638 ± 0.00066 


0.000 0.000 0.000 0.000 1.000 
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The precise measurements of LEP and SLC and their correlations |25 cii'g summarized 
in Table IHl Most of these Z observables would be modified in case of an additional New 
Physics contribution to hadronic Z decays. Let e^p be the ratio of this new partial width 
Pnp to the total decay width of the Z without this new contribution. As was shown in 
Ref. 2 total width Pz, the ratio Ri of the hadronic to the leptonic branching 

fractions, and the peak cross section a^^^ are modified as follows, 



Tz 

Re 



Pz (i + , 

Re (l + , 



[rz + P 



NP 



a 



had 



(l - 0.57£ 



had 
NP 



[(Fhad + Tnp) /Pf] 
12n Pee(rhad + Tnp) 

m| (Pz + Pnp)2 



(3) 
(4) 

(5) 



In Ref. the new hadronic decay channel considered was flavour-democratic. The 
individual branching fractions into the different quark flavours were therefore not modified 
by this new contribution. In the case of a sbottom pair production with hadronic R-parity- 
violating decays into light quarks exclusively, the ratio of the bb branching ratio to the 
hadronic branching ratio, i?b, is also modified according to 



Rh 



Rh 



lA3e 



had 
NP 



i^hh/ (rhad + Tnp)] 



(6) 



while (gv/gA) remains untouched. 

These observables would also be modified by the virtual corrections arising from the 
New Physics responsible for the additional hadronic contribution. As in Ref. |2Z|; the 
value of e^p was fitted to the measurement of the five observables together with the 
generic contribution of these virtual effects. The result is 



_had 
-NP 



(-0.56 ±0.80) X 10- 



(7) 



which corresponds to an additional hadronic contribution of 



<d("^z) = -24±36pb. (8) 

It allows a 95% CL. upper limit of 56 pb to be set on the cross section, at the Z peak, of 
any additional hadronic contribution to the Z decays into light quarks only. The resonant 
contribution of the sbottom pair production cross section with = 6 GeV/ is 
shown in Fig. El as a function of the mixing angle cos6'mix between the two sbottom states 
bi and Br, superpartners of the left-handed and right-handed bottom quarks, respectively. 
For cos6'mix — 0.39, the coupling between the Z and the lighter sbottom vanishes. 

For = 6 GeV/c^, the Z data allow all values of cos^mix below 0.22 and above 0.52 
to be excluded at the 95% confidence level. These data are therefore incompatible with a 
light sbottom pair production, unless the coupling to the Z is negligibly small. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos© 

mix 

Figure 2: The resonant contribution of the sbottom pair production cross section with — 6 GeV/c^, 
at a/s — mz, as a function of cos^mix (full curve). The dash-dotted line indicates the 95%C.L. upper 
limit on this cross section when the sbottom decays into light quarks exclusively. 



2.3 The LEP 2 data 

The preliminary LEP 2 hadronic cross section data were taken from Ref. [21] • The mea- 
sured cross sections cxhad and the standard model predictions cxth are summarized in Ta- 
ble [7| These data are displayed in Fig. El and the correlation matrix is given in Table |H1 
When no systematic uncertainties are assigned to the theoretical prediction, the average 

Table 7: The hadronic cross section, dhad, measured at the twelve LEP 2 centre-of-mass energies, and 
the predictions in the standard model, uth- These data are still preliminary. 



(GeV) 


o-had (pb) 


o"th (pb) 


130 


82.1 ±2.2 


82.8 


136 


66.7 ±2.0 


66.6 


161 


37.0 ± 1.1 


35.2 


172 


29.23 ±0.99 


28.74 


183 


24.59 ± 0.42 


24.20 


189 


22.47 ±0.24 


22.16 


192 


22.05 ±0.53 


21.24 


196 


20.53 ±0.34 


20.13 


200 


19.25 ±0.32 


19.09 


202 


19.07 ± 0.44 


18.57 


205 


18.17 ±0.31 


17.81 


207 


17.49 ±0.26 


17.42 



ratio appears to exceed the prediction by (1.5 ± 0.9)%, i.e., by 1.7 standard deviations. 
This excess, although not significant, is compatible with, and actually slightly larger 
than an additional light sbottom pair production with cos^mix = 0.39. The latter would 
amount to about 1% of the total cross section in this centre-of-mass energy range. 
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Figure 3: Ratio (a) and difference (b) of the hadronic cross section measurements and the standard 
model prediction as a function of the centre-of-mass energy, for the LEP2 data. The dash-dotted line 
indicates the standard model prediction, and the dashed curve the additional contribution of sbottom 
pair production with mj^ = 6 GeV/c^ and cos^mix = 0.39. 



Table 8: The correlations between the twelve LEP 2 centre-of-mass energy bins (v^ in GeV). 
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3 Global fit 



When no systematic uncertainties are assigned to the standard model prediction, the data 
can be combined in a global negative log-likelihood £(cos6'mix, «) as follows, 

1 ^ 

/:(cos6'mix, ") = 2 H ^i^ij^^j wi^h Ai = (Thad,i - Wth,i + acrNp,i(m^, cos6'mix)] , (9) 

where S is the covariance matrix of the N (= 28) measurements of PEP, PETRA, TRIS- 
TAN, LEPl, SLC and LEP2 as compiled in Section |2l 6'mix is the mixing angle in the 
sbottom sector and a is an arbitrary normalization constant of the sbottom pair produc- 
tion cross section, cjNp.i- The likelihood is then minimized with respect to cos^mix and to 
a to find the best fit to the data. A fitted value of a compatible with unity and incom- 
patible with 0. would be the sign of New Physics, while a value compatible with 0., but 
incompatible with 1., would allow this New Physics to be excluded with a certain level of 
confidence. (This same technique can be applied for any kind of New Physics leading to 
hadronic final states in e^e~ collisions.) 

For a = 1 and m^^ = 6 GeV/c^, the negative log-likelihood is displayed in Fig. 0^ as a 
function of cos^mix. Not surprisingly, the Z peak data f Section 12. 2|) constrain the coupling 
of the sbottom to the Z to be vanishingly small, cos6'inix = 0.39 ± 0.07. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 




-1 -0.5 0.5 1 1.5 2 2.5 3 

a 

Figure 4: The negative log-likelihood with m^^ — 6 GeV/c^ (a) as a function of cos^mix for a = 1; and 
(b) as a function of a with cos6'mix = 0.39 for the combined data (full curve), PEP/PETRA (dashed 
curve), TRISTAN (dotted curve) and LEP 2 (dash-dotted curve). 
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The value of the mixing angle was therefore fixed to cos^mix = 0.39. The combined 
negative log-likelihood and those for PEP/PETRA, TRISTAN and LEP 2 data are shown 
in Fig. as a function of a. (For LEP 1 and SLC, the likelihood does not depend on 
a, because of the vanishing sbottom cross section for cos6'inix = 0.39.) The values of 
a for which the different negative log-likelihood functions are minimized are indicated 
in Table IHl together with the corresponding 68% confidence intervals and the 95% CL. 
upper limits on a. (This one-sided upper limit is the a value for which the negative 
log- likelihood increases by 1.64^/2 with respect to the minimum.) 

Table 9: The values a^^^ for which the negative log-hkehhood is minimized in PEP/PETRA, TRISTAN 
and LEP 2 data, and in the combination, together with the 68% confidence intervals and the 95%C.L. 
upper limits, ags, for cos6'mix — 0.39 and m^^ = 6GeV/c^. 



Data 




«95 


PEP/PETRA 


0.45 ±0.30 


0.94 


TRISTAN 


0.21 ±0.39 


0.85 


LEP 2 


1.32 ±0.74 


2.52 


All 


0.45 ±0.23 


0.82 



As was already alluded to in Section 12.11 the lower energy data do not favour the 
sbottom hypothesis [a = 1). They are, instead, compatible with the standard model 
{a = 0) within one standard deviation or thereabout. A slight excess in the LEP 2 data, 
at the 1.7a level f Section 12. 3|) . translates as such to the combined result. The latter, 
however, excludes the sbottom hypothesis with = 6GeV/c^ at more than 95%C.L., 
when no systematic uncertainty is assigned to the standard model prediction. 

The main sources of uncertainty for the theoretical prediction of the e"'"e~ —>■ qq cross 
section are (i) the knowledge and the running of the strong coupling constant as', (ii) 
the running of the electromagnetic coupling constant oqed; and (in) the theoretical ac- 
curacy of the prediction from the ZFITTER program. As in Ref. 27j, the values and the 
uncertainties of the strong and electromagnetic coupling constants were taken to be 

a5(mz) = 0.1183 ±0.0020 [2i and a(mz)"^ = 128.95 ± 0.05 [HI, (10) 

leading to uncertainties in the hadronic cross section prediction of 0.15% and 0.08%, 
respectively. The missing higher orders in ZFITTER are estimated to contribute another 
0.1%. These numbers add quadratically to a total systematic uncertainty r^th of the order 
of 0.2%, in agreement with the estimate of Ref. [21] (?7th = 0.26%) for LEP 2 data. 

If this common systematic uncertainty is assumed to have a Gaussian probability 
density function, the negative log-likelihood can be modified as follows, to account for the 
full correlation between all centre-of-mass energies: 

^ = o E KSij'^'j + ^ with = (Thad,^ - [(1 ± Pth)ath,i ± a^NP.^] , (11) 

^ i,j=l '^^th 

where pth is the actual theoritical bias of the standard model prediction, to be fitted from 
the data. 
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It is reasonable, however, to take into account the non-Gaussian nature of uncertainties 
of theoretical origin. For example, the missing higher orders in ZFITTER may turn into 
a bias of —0.1%, 0% or 0.1% with an equal probability. (In fact, the least likely value is 
certainly 0%, as missing orders are expected to contribute a finite amount to the cross 
section.) Similarly, the uncertainty on the absolute value of as'(mz) is dominated by 
theory, and cannot be considered as Gaussian. It is therefore probably more adequate to 
assume a probability density function as displayed in Fig. i.e., flat between — ?7th and 
+rjth^ and with a Gaussian shape outside this interval. 



1.5 - 



1 - 



0.5 - 




-1 -0.5 0.5 1 

Relative bias (%) 

Figure 5: Probability density function for the conventional (Gaussian) systematic uncertainty treatment 
(dashed curve) and suggested here instead (full curve) to account for the non-Gaussian nature of theory 
uncertainties, with 77th = 0.2%. 

The likelihood was therefore further modified by changing the Pth/2^th term to 

(pth + r/th)V2^th if Pth < -r?th, (12) 

if -r^th < Pth < ^?th, (13) 
(pth -^/th)V2^th if Pth > r?th. (14) 



This negative log-likelihood was then minimized with respect to the theoretical bias 
Pth, for each value of a, with Gaussian and non-Gaussian uncertainties. The result is 
displayed in Fig. El in the two configurations as a function of a, for cos6'mix = 0.39 and 
= 6GeV/c^. The values of a for which the negative log-likelihood is minimized are 
indicated in Table ^1 together with the corresponding 68% confidence intervals and the 
95%C.L. upper limits on a. 
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Table 10: The values amin for which the combined negative log-hkeUhood is minimized for Gaussian and 
non-Gaussian uncertainties, together with the 68% confidence intervals and the 95%C.L. upper limits, 
a95, for cos6»mix = 0.39 and = GGeV/c^. The fit results for PEP/PETRA, TRISTAN and LEP2 
(apETRA, QfTRiSTAN and aLEP2) are also given. 



Uncertainties 




"95 


ttPETRA 


^TRISTAN 


aLEP2 


Gaussian 
Non- Gaussian 


0.45 ± 0.24 


0.85 
0.92 


0.45 ±0.35 
0.59t°J^ 


0.16 ±0.47 

U.UO_o.48 


1.68 ± 1.02 
1.87 ± 1.02 



It can be seen that the upper hmit on a depends very httle on the way the common 
systematic uncertainties are dealt with. The most conservative approach is chosen here 
to derive the final results. 
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4 Results 



The same procedure was repeated by varying the sbottom mass from to 12GeV/c^. 
For each mass, the 95% C.L. upper hmit on a was determined as explained above. A 
sbottom with a given mass is excluded if this upper limit is smaller than unity. Figure [3 
shows the 95% C.L. upper limit on a for cos^mix = 0.39 as a function of the sbottom 
mass, with Gaussian and non-Gaussian uncertainties. (In the latter configuration, the 
non-Gaussian nature of the likelihood was taken into account in the determination of the 
limit.) Sbottom masses below 7.5GeV/c^ are excluded at the 95% confidence level. 



1.4 



a 
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u 
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u 
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1.2 - 
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Non-Gaussian uncertainties 
Gaussian uncertainties 



95%C.L. Excluded 



Ref. V 
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..:■////////// 





0.8 




COS© 



:0.39 



0.6 



10 



12 



Sbottom mass (GeV/cO 



Figure 7: The 95% C.L. upper limit on a as a function of the sbottom mass, with non-Gaussian (full 
curve) and Gaussian (dashed curve) common systematic uncertainties. Also shown are the predictions of 
the model of Ref. 0], now excluded by this analysis. 



Because cos^mix is very much constrained by the Z peak data, the upper limit on a 
is expected to be smaller than that shown in Fig. [7| for any other value of the mixing 
angle. As a check, the procedure was repeated again by varying cos^mix from to 1, with 
non-Gaussian uncertainties. The resulting sbottom mass lower limit is shown in Fig. |H1 
as a function of cos6'mix, and is indeed at least 7.5GeV/c^ over the whole range. (The 
region excluded by LEP 2 data at large values of cos ^mix is probably over-optimistic, as 
four-jet events - expected from such heavy sbottom pair as well as W pair production - 
are rejected from the qq event samples selected above the WW threshold.) 

It is worth mentioning that the presence of a light sbottom would slow down the 
running of 05 with the centre-of-mass energy. (It would be even more so with an additional 



16 



80 



o 



o 



U 

ON 



Hadronically decaying sbottom 

Exeliitsion by PEP, PETRA, TRISTAN i 
Exclusion by the sole LEPl data 




Figure 8: Absolute 95%C.L. lower limit on m^^ as a function of cos^mix, for hadronically decaying 
sbottoms. The hatched area is excluded at 95%C.L. The dashed line shows the exclusion achieved with 
the sole Z peak data. 



light gluino.) Starting from the value accurately measured in r decays [30j, (the only 
measurement not affected by a sbottom heavier than 2 GeV/c^ and lighter than 5.5 GeV/c^, 
and corresponding to as{mz) = 0.121 ±0.003 in the standard model), this slower running 
would lead to values of 0:5 larger than assumed in this letter, at all centre-of-mass energies. 
The total New Physics contribution (from the direct sbottom production and the increase 
of as) would further increase the effect on the total hadronic cross section expected at 
PEP, PETRA, TRISTAN, SLC and LEP. The 7.5 GeV/c^ lower limit on the sbottom mass 
is therefore probably very conservative. 



5 Conclusion 

The e+e~ —>■ hadron cross section data collected well above the bb resonances have been 
compiled and analysed to search for an anomalous production of hadronic events. Alto- 
gether, the PEP, PETRA, TRISTAN, LEP 1, SLC and LEP 2 data allow a light sbottom 
decaying hadronically to be excluded at 95% CL. for any mixing angle, if its mass is below 
7.5 GeV/c^. When combined with the result of Ref. |Sj in which a stable sbottom with 
mass below 92 GeV/ is excluded, this analysis definitely invalidates the model of Ref. ^ 
with a 12-16 GeV/ gluino and a 2-5.5 GeV/ sbottom. 
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Note added 

This work has been primarily motivated by the "apparent excess" reported in Ref. With 
the coUaborative help of the author, the excess was found to be an artifact of duphcated, missing 
and over-corrected data in the computer-readable files of the Review of Particle Physics [221 123j 
augmented by an incorrect interpretation of the Z peak data in Ref. 31 . The aforementioned 
computer-readable files are being updated to include the work described in this letter. 
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